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Composición y estructura de asociaciones de insectos 
acuáticos de arroyos de selva nublada de Yungas del NO 
Argentino 


a RESUMEN. Treinta y tres ambientes lóticos en las selvas nubladas 
de топтаћа de las Yungas del NO Argentino, fueron muestreados tanto 
en ambientes no modificados como alterados por actividades humanas. 
Insectos acuáticos de 143 taxones en 55 familias fueron colectados. El 
análisis de agrupamientos sugirió que la altura es una de las variables 
principales en la estructuración de las comunidades de insectos en estos 
arroyos, y la importancia de esta variable fue confirmada mediante un 
ordenamiento no-métrico multi-dimensional (NMS); los parámetros 
ambientales que mejor se correlacionaron con la ordenación fueron: 
altura, temperatura del agua, latitud y variables del canal (ancho, 
porcentaje de detritos leñosos grandes y pequeños, de bancos excavados, 
de piedras y grava gruesa). Procedimientos de permutación de respuestas 
múltiples (MRPP), mostraron que los arroyos en áreas bien conservadas 
difieren significativamente en su composición de los arroyos en áreas 
modificadas. La proporción de individuos de Elmidae y Plecoptera, y el 
número de taxones de Trichoptera, fueron los métricos biológicos mejor 
correlacionados con el gradiente de alteración ambiental local, sugiriendo 
que un índice ‘EIPT’ podría ser un componente útil para la evaluación 
del estado ecológico de estos ambientes. Los análisis de indicadores de 
especies, identificaron algunos indicadores potenciales de la condición 
de los arroyos y de los factores de alteración que los afectan. 


PALABRAS CLAVE. Diversidad. Estructura de la comunidad. Especies 
indicadoras. Índices biológicos. Elmidae. 


Ш ABSTRACT. Thirty three lotic environments in the Yungas 
mountain cloud forest of NW Argentina were sampled both in undisturbed 
forest areas and sites altered by human activities. Aquatic insects of 143 
taxa in 55 families were collected. Cluster analysis suggested altitude as 
one of the main structuring variables of aquatic insect communities in 
these streams, and its importance was confirmed by non-metric multi- 
dimensional scaling (NMS); the environmental parameters measured that 
were best correlated with the ordination were altitude, water temperature, 
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latitude and channel variables (width, percentage of large and small 
woody debris, of undercut banks, cobble and coarse gravel). Multi 
response permutation procedures (MRPP) showed streams in well 
preserved areas to significantly differ in their composition from streams 
in disturbed areas. Proportion of Elmidae and of Plecoptera individuals 
and number of Trichoptera taxa were the biological metrics best correlated 
with the local disturbance gradient, suggesting that an ‘EIPT’ index could 
be a useful component in the evaluation of the ecological status of these 
environments. Indicator species analyses identified some potential 
indicators of stream condition and disturbance factors affecting these 


streams. 
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INTRODUCTION 


The Yungas are cloud forests extending 
from the south of Venezuela into NW 
Argentina along the Andean mountain range. 
In Argentina they are distributed 
discontinuously within the provinces of Salta, 
Jujuy, Tucumán and Catamarca, along part 
of the Subandean chains, and they represent 
one of the most species-rich biogeographic 
provinces (Brown, 1995). Climate is warm 
and humid to sub-humid, and altitude 
ascends Кот 300 to 2400 m.a.s.l. 
Temperature and humidity vary in relation to 
altitude, latitude and aspect of slopes. 
Average yearly precipitation ranges from 900- 
1000 mm, reaching 1300 mm in some areas. 
Rains are concentrated within 5 to 6 months 
in summer, and during colder months, fog 
partly compensates for lack of rainfall (Burkart 
et al., 1994). With increasing altitude, 
different altitudinal floors or vegetation levels 
are recognizable: warm and humid foothill 
jungle; temperate-warm and humid mountain 
jungle; and temperate and humid mountain 
forest with frequent winter frost. 

The Yungas have become intensively 
fragmented during the last few decades due 
to ongoing agricultural practices (Reboratti, 
1989). Human activities such as conversion 
of forest to agricultural or residential areas, 
logging, mining, road construction and 
recreation, modify reach and channel 
variables which are well known factors 
structuring aquatic insect assemblages (Allan, 
2004). Sponseller et a/. (2001) found that 
macroinvertebrate indices were more closely 


related to land-cover patterns evaluated at a 
local 200 m corridor scale, suggesting that 
streamside development effectively alters 
assemblage structure. In Patagonia, loss of 
riparian cover as a consequence of 
overgrazing and other land use practices was 
shown to have negative effects on aquatic 
environments (Miserendino, 2004). Several 
studies on benthic macroinvertebrates of lotic 
environments have been conducted in the 
Yungas of NW Argentina during the last years 
(Fernández & Fernández, 1995; Domínguez 
& Fernández, 1998; Romero & Fernández, 
2001; Fernández et al., 2001; 2002; 2006; 
Fernández & Molineri, 2006), reflecting the 
rising concern for loss of biodiversity and 
alteration of landscape in the region. In order 
to evaluate the ecological status (Prat et a/., 
1999) of Yungas streams, studies integrating 
information from their riparian areas into the 
analysis are needed (Fernández et a/., 2006). 

Sites of present study are within the 
Bermejo basin in Salta and Jujuy provinces, 
where lotic environments (streams and small 
rivers) included in the three altitudinal floors 
of the Yungas were sampled, both in 
protected and non-protected areas. My goals 
were to document composition and structure 
of insect assemblages of Yungas streams in 
relation to natural gradients, evaluate the 
response of stream condition to a disturbance 
gradient at the local scale, and explore the 
potential for indicators of stream condition 
and disturbance factors. | dedicate this paper 
with affection to Dr. A.O. Bachmann, our 
father of aquatic entomology, who has always 
been available to provide advice on 
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biological, systematic or nomenclatorial 
questions when answers from books were not 
forthcoming. 


MATERIAL AND METHODS 
Study area and sampling protocol 


Thirty three lotic environments located in 
both protected and non-protected areas (Fig. 
1) within four of the Yungas forest fragments 
in Salta and Jujuy provinces were sampled 
between August 2005 and May 2006. In each 
locality three samples each lasting 15 minutes 
were taken and later pooled for 
standardization, encompassing all 
microhabitats of different depths, substrate 
and current along three transects separated 
by 10 m across the stream. Aquatic insects 
were collected with aquatic nets of 20 cm 
diameter and mesh size of 300 um, and 
manually by visually inspecting rocks, 
substrate and aquatic vegetation. Specimens 
are preserved in the collections of the Museo 
de Ciencias Naturales de Salta and Instituto 
y Fundación Miguel Lillo, Tucumán 
(Argentina). Biotic and abiotic data were 
recorded at each station (Appendix 1) 
including qualitative (visual estimates) and 
quantitative in-stream, channel and riparian 
characteristics, as well as evidence of habitat 
disturbance (presence of roads, garbage, 
cattle, crops, logging, channelization) at the 
local reach scale, extending 100 m in width 
on each bank and 300 m in length upstream. 
Stations were classified (Table |) according 
to altitudinal level as follows: 1. foothill 
rainforest (ca. 300-900 m), 2. mountain 
rainforest (ca. 900-1500 m) and 3. mountain 
forest (ca. 1500-2400 m), as well as according 
to the disturbance level observed at the local 
reach scale (riparian area and channel) as 
follows: 1. well preserved: vegetative cover 
of tall trees (about 25 m high), thick riparian 
vegetation with по evidence of habitat 
alteration, and disturbed, in three categories: 
2. preserved: trees 15-20 m high, bushes 
comprising large part of the cover, as indirect 
evidence of habitat alteration; 3. modified: 
vegetative cover with bushes making up for 
most of it, direct evidence of habitat 


alteration; 4. highly modified: vegetative 
cover scarce and riparian forest deteriorated, 
direct evidence of habitat alteration. 


Data analysis 
Characterization of the community 


Taxonomic identifications were made to 
genus and species level when possible. In 
cases of doubtful identification, taxa were 
either combined (ѓе. Massartellopsis/ 
Meridialaris in Leptophlebiidae and 
Banyallarga/Phylloicus in Calamoceratidae) 
or identified to family (ѓе. Chironomidae in 
Diptera). Information was arranged in a 
matrix of 143 taxa by 33 sampling stations 
(Appendix 2) based on which biological 
metrics of number of individuals and of taxa 
(at taxonomic level from Appendix 2) per 
order, and proportions of both to the total 
were calculated for each station. 

Three diversity indices were estimated 
based on Whittaker (1972): alpha diversity, 
calculated as the average specific richness per 
locality; beta diversity, a measurement of the 
heterogeneity of the data, calculated as the 
ratio between total number of species and 
average number of species; and gamma 
diversity, or diversity at landscape level, 
calculated as total number of species across 
all localities. Specific richness, evenness and 
diversity (Shannon and Simpson indices) were 
also calculated for each site. 


Natural factors affecting the community 


In order to avoid the masking effect that 
pollution could have on the analyses of 
natural factors affecting community 
composition, potentially polluted sites were 
excluded from these analyses. Polluted sites 
were identified based on the BMWP biotic 
index adapted for the region by Domínguez 
& Fernández (1998). Since this index is based 
on presence/ absence of families of benthic 
macroinvertebrates, non-benthic insects were 
excluded from the calculation, which 
included only benthic insects and Crustacea 
(no Oligochaeta, Mollusca or Hydracarina 
were collected). Although BMWP values 
obtained might be underestimated (for 
exclusion of some non insect taxa) and not 
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Fig. 1. Study area. Hatched: Yungas forest; Grey polygons: Protected areas; 1: National Park Baritú; 2: 
Private Reserve El Pantanoso; 3: National Park Calilegua; 4: National Park El Rey; Contour: Biosphere 
Yungas Protected area. Coded squares: Localities: Yu1: stream 1 affluent to Baritú River; Yu2: Baritü 
River; Yu3: stream 2 affluent to Baritú River; Yu4: stream 1 crossing road Baritú-Lipeo; Yu5: stream 2 
crossing road Baritú-Lipeo; Yu6: stream affluent to Cayotal River; Yu7: Cayotal River at Las Juntas; Yu8: 
stream affluent to Lipeo River; Yu9: stream 1 crossing road Lipeo-Los Toldos; Yu10: stream 2 crossing 
road Lipeo-Los Toldos; Yu11: stream 3 crossing road Lipeo-Los Toldos; Yu1 2: effluent river from Dique 
Itiyuro; Yu13: Pantanoso River; Yu14: Zora River; Yu15: San Francisco River; Yu16: stream at Route 6 
to Palma Sola; Yu17: Tres Cruces Stream; Yu18: stream in Lesser; Yu19: La Sala Stream 1; Yu20: La Sala 
Stream 2; Yu21: stream on trail to Pozo Verde; Yu22: stream in La Angostura; Yu23: stream in Quebrada 
de San Lorenzo; Yu24: Vaqueros River; Yu25: stream in Las Costas; Yu26: Las Nieves Stream; Yu27: 
Yuto River; Yu28: stream in Quebrada de Tilián; Yu29: Zanjón Seco Stream; Yu30: stream 20 km SE to 
Isla de Cañas; Yu31: stream 15 km SE to Isla de Cañas; Yu32: stream 5 km SE to Isla de Cañas; Yu33: 
Anta Muerta Stream. 
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comparable to values of other studies, they 
do allow the identification of non- polluted 
sites. 

Composition of aquatic insect 
communities was analyzed using multivariate 
cluster and ordination techniques with the 
program PC-ORD (McCune & Grace, 2002). 
Both multivariate analyses were based on 
abundance data (values transformed as log 
(x+1) to homogenize the variance, reduce the 
effect of total quantity emphasizing relative 
quantities, and equalize relative importance 
of common and rare species) from the 30 
non- polluted sites (water quality class 1, 
Table |). 

For the cluster analysis, Sorensen (Bray- 
Curtis) distance coefficient was used, with 
flexible beta as linkage method at a value of 
B = -0.25. The resulting dendrogram was 
based on Wishart’s objective function 
converted to a percentage of remaining 
information. A value of 25% of remaining 
information was selected to prune the 
dendrogram. 

Because of its ability to extract information 
from nonlinear relationships, non-metric 
multidimensional scaling (NMS; Mather, 
1976) was chosen as ordination method. Its 
use and interpretation differ from other 
ordination methods in that: 1. for a given 
number of dimensions the solution for a 
particular axis is unique; 2. the appropriate 
number of dimensions (axes) is determined 
graphing or tabulating final stress versus 
number of dimensions, choosing a number 
of axes after which reduction in stress is small; 
and 3. axis numbers are arbitrary, so that the 
percentage of variance on a given axis does 
not necessarily form a descending series with 
increasing axis number. Sorensen was chosen 
as distance coefficient. Forty runs were 
carried out with real data and 50 with random 
data (Monte Carlo test) starting from a random 
configuration, and with a possible maximum 
of 6 axes and 400 iterations. Final instability 
was Calculated as standard deviation in stress 
over the preceding 15 iterations (value of < 
10% indicates a stable solution; McCune & 
Grace, 2002). Proportion of variance 
represented by ordination axes was 
calculated by the correlation (determination 
coefficient 1?) between Euclidean distances 


in the ordination space and distances in the 
original space. Distances in the original space 
were calculated with the same distance 
measure used in the NMS analysis. Since pH 
data were not available for all localities, the 
NMS analysis was run twice, once including 
all localities and the second time excluding 
localities lacking pH data. 

Quantitative environmental variables 
(transformed as log (х+1) to make the units 
of attributes measured in different scales 
comparable) listed in Appendix 1 were 
correlated with NMS ordination axes and 
overlaid onto the ordination diagrams as joint 
plots (diagrams of radiating lines whose angle 
and length indicate direction and strength of 
relationships between variables and 
ordination axis) where only variables from 
the environmental matrix with an 2 > 0.20 
were represented. Groups obtained with the 
cluster analysis were also overlaid onto the 
NMS ordination to aid in the interpretation 
of their relationships. 

An indicator species analysis based on the 
method described by Dufréne & Legendre 
(1997) was conducted to identify potential 
indicators for the three altitudinal floors of 
the Yungas. This method combines 
information about species abundance and 
frequency of occurrence in each group. A 
perfect indicator for a particular group 
(indicator value of 100) must be faithful 
(always present) and exclusive to the group 
(never occurring in other groups). Statistic 
significance of the indicator values was 
established with a Monte Carlo test (with 
1000 permutations). 


Analysis of disturbance 


Spearman rank correlations between 
disturbance level gradient (from well 
preserved, class 1, to highly modified, class 
4) and natural environmental variables were 
performed in order to determine if they 
covary (values transformed as log (х+1)). 

Sites in well preserved areas according to 
the local disturbance level (11 localities, 
disturbance level 1 in Table |) were compared 
to sites in disturbed areas (22 localities, 
disturbance levels 2-4 in Table I), using the 
following methods: 
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Fig. 2. Cluster analysis dendrogram based on similarity of aquatic insect composition among sampling 
stations. Numbers indicate groups formed when pruning the dendrogram at a level of 25 % of remaining 
information (vertical line). Altitudinal floor is superimposed to show correlation with stations. 


- A MRPP (multiple response permutation 
procedure) test. This method provides a 
multivariate non-parametric test of 
differences among two or more groups based 
on the analysis of a distance matrix. Sorensen 
was used as distance coefficient. Delta 
(average weighted distance within a group; 
lower delta value indicates better cohesion 
within a group) was calculated according to 
the procedure detailed by Mielke & Berry 
(2001). The statistic of this method, T= 
observed - expected delta/ М variance of delta, 
describes separation between groups (the 
more negative the value of 7, the larger the 
separation between groups), А = describes 
homogeneity within each group compared 
to one due to chance, and pz the probability 
of obtaining a delta as high as or higher than 
observed given the distribution of possible 
deltas (probability that observed difference 
is due to chance). 

- Comparison of richness, diversity 
indexes and biological metrics through 
Student t tests, illustrated with side by side 
boxplots. 

- Spearman Rank correlations between the 
disturbance level gradient and biological 


indexes and metrics (values transformed as 
log (x+1)). 

- Indicator species analyses based on 
Dufréne & Legendre's (1997) method 
(described in previous section) to identify 
potential indicators for sites in well preserved 
and disturbed areas, and for groups of sites 
defined by presence or absence of different 
factors of habitat disturbance (cattle, crops, 
logging, etc., listed in Appendix 1). 


RESULTS 
Characterization of the community 


A total of 143 aquatic insect taxa was 
found (Appendix 2). Alpha diversity was 20.2, 
beta diversity 7.07, and gamma diversity 143. 
Species richness varied from 7 to 41 across 
localities, and localities with the highest 
richness values corresponded to well 
preserved environments within protected 
areas (Yu19, Yu18 and Yu20 in National Park 
El Rey, and Yu13 in El Pantanoso Reserve). 
Localities with the highest diversity values 
were Yu19 and Yu21 (National Park El Rey), 
Yu 13 (El Pantanoso Reserve), Yu3 and Yu8 
(National Park Baritú), and Yu31 (Arroyo 15 
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Fig. 3. Joint plots showing relationships between environmental variables and axes of NMS ordination 
of localities. 3 D solution reached after 122 iterations, with a final stress of 15.07 and final instability of 
105. 3.1: Axes 1 and 2; 3.2: Axes 1 and 3; 3.3: Axes 2 and 3. 
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Km E to Isla de Cañas, not protected) (Table 1). 
Natural factors affecting the community 

Thirty localities were classified as ‘non- 
polluted’ according to the BMWP index 
approximation (water quality class 1, Table 
|) and were included in the analyses 
evaluating influence of natural factors. 

Hierarchical cluster analysis identified five 
groups with 25 96 of remaining information 
among members for each group. The cluster 
including groups 1, 2 and 3 joined all 
localities encompassed within the mountain 
forest altitudinal floor, while the cluster 
including groups 4-5 encompassed those 
within the foothill rainforest altitudinal floor 
(Fig. 2), suggesting that altitudinal gradient is 
an important factor in shaping aquatic insect 
assemblages in the Yungas. 

A three dimensional NMS final solution 
was reached after 122 iterations, with a final 
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stress of 15.07, final instability of 10°, and a 
proportion of randomized runs with stress 
lower than or equal to observed stress of 
0.0196. Variance represented by the three 
ordination axes was of 30.7, 18.3 and 30.5 
(cumulative variance of 79.6) respectively. 
Correlation of environmental variables 
with ordination axes (Table II) shows that the 
most important environmental variables 
explaining community composition include 
channel variables of percentage of coarse 
gravel in the substrate, and large woody 
debris and undercut banks in the channel 
along ordination axis 1 (Figs. 3.1, 3.2), while 
altitude and wet channel width, percentage 
of cobble in the substrate and of small woody 
debris in the channel significantly contributed 
to axis 2 (Figs. 3.1, 3.3), and altitude, latitude, 
and water temperature to axis 3 (Figs. 3.2- 
3.3). Biological groups obtained from cluster 





Richness 
p 0.0084 


% Ephemeroptera Indv 


Fig. 4. Side by side box plots of Shannon diversity and some biotic metrics comparing median, quartiles 
(25th and 75th percentiles) and 9596 confidence intervals of sites in well preserved (WP, М=11) and 
disturbed (Di, N222) areas. p values obtained with Student ttests. 
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Table I. Locality classification and summary statistics. BMWP (Domínguez & Fernández, 1998: >40= 
non-polluted (water quality class 1); 30-40= slightly polluted (water quality class 2); 20-30: polluted 
(water quality class 3). See text for explanation of other codes. 























































































































Locali Altitudinal | Disturbance BMWP Water Quality Number of ; Shannon Simpson 
ocality fl level irid а А Richness | Evenness Diversit Diversit 
oor evel ndex ass specimens iversity iversity 
Yul 2 72 1 58 18 0.89 2.58 0.90 
Yu2 2 2 119 1 120 27 0.85 2.80 0.91 
Yu3 2 89 1 118 26 0.90 2.94 0.93 
Yu4 3 55 1 62 16 0.88 2.44 0.88 
Yu5 2 2 38 2 30 7 0.78 1.52 0.72 
Yu6 2 3 35 2 48 11 0.88 212 0.85 
Yu7 3 2 62 1 66 17 0.89 2.53 0.90 
Yu8 3 55 1 49 17 0.94 2.66 0.92 
Yu9 3 3 43 1 34 10 0.77 1.77 0.74 
Yu10 1 69 1 59 20 0.87 2.60 0.89 
Yul 1 70 1 78 2 0.80 2.42 0.85 
Yu12 1 4 29 3 70 9 0.86 1.88 0.81 
Yu13 1 104 1 199 37 0.83 3.00 0.92 
Yul4 2 3 79 1 331 2 0.47 1.42 0.53 
Yu15 2 3 40 1 175 16 0.72 1.99 0.82 
Yul6 2 2 59 1 241 20 0.53 1.59 0.63 
Yu17 2 40 1 131 23 0.74 2.31 0.80 
Yu18 2 3 70 1 127 3 0.84 2.89 0.92 
Yu19 2 108 1 225 4 0.87 3.21 0.94 
Yu20 3 94 1 245 3 0.72 2.47 0.81 
Yu21 2 92 1 57 26 0.89 2.91 0.93 
Yu22 3 3 45 1 54 15 0.87 2.35 0.86 
Yu23 3 2 44 1 93 15 0.72 1.96 0.79 
Yu24 2 3 80 1 135 20 0.79 2.37 0.86 
Yu25 2 3 48 1 35 10 0.74 1.71 0.73 
Yu26 1 2 69 1 80 22 0.78 2.41 0.86 
Yu27 1 2 44 1 75 17 0.80 2.26 0.85 
Yu28 1 2 29 1 127 22 0.79 2.45 0.86 
Yu29 1 3 95 1 57 12 0.80 1.99 0.81 
Yu30 1 3 17 1 119 19 0.80 2.37 0.85 
Yu31 2 3 42 1 97 28 0.88 2.95 0.93 
Yu32 2 4 92 1 71 20 0.87 2.62 0.90 
Yu33 2 4 61 1 90 22 0.71 2.21 0.78 
Mean 63.65 104.54 19.69 0.80 2.35 0.84 
SD 25.54 70.63 7.87 0.10 0.46 0.09 



































Table Il. Correlations of environmental variables with NMS ordination Axes. Variables with an 7? 


larger than 0.20 shown in joint plots of Fig. 3 are highlighted in bold. г Pearson Product-Moment 
coefficient; tau: Kendall’s tau rank coefficient. 





Variables tau 

Altitude 0.671 | «.000 : -0.402 
Latitude -0.008 | 0.968 0.402 
Water temperature -0.34 | 0.066 0.408 
96 cobble 0.516 | 0.004 -0.02 
% coarse gravel 0.34 | 0.066 -0.248 
% mud -0.445 | 0.014 0.027 
Width -0.49 | 0.004 -0.067 
% filamentous algae -0.004 0.425 0.352 
% macrophytes -0.31 | -0.165 | 0.381 0.237 
95 Ig. woody debris 0.568 | 0.001 | 0.323 0.484 0.252 | 0.179 | 0.063 0.184 0.158 | 0.403 | 0.025 0.142 
% sm. woody debris | 0.095 | 0.614 | 0.009 | 0.123 0.468 | 0.009 | 0.219 | 0.368 | -0.121 | 0.527 | 0.015 | -0.068 
96 undercut banks 0.517 | 0.003 | 0.267 | 0.412 | -0.221 | 0.240 | 0.049 -0.2 0.079 | 0.678 | 0.006 | 0.055 
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Table III. Indicator taxa for altitudinal floors obtained by combining abundance and frequency data 
according to Dufréne & Legendre's (1997) method. Only taxa with a statistically significant indicator 
value (larger than expected by chance, p « 0.05) are shown. Indicator values range from 0 to 100 
(perfect indicator). 
































Taxa Altitudinal floor Indicator value Mean SD p 

Americabaetis 1 75.7 37.3 11.12 0.002 
Rhagovelia (Neorhagovelia) 1 65 32.3, 9.95 0.012 
Chironomidae 1 54.4 34.5 7.67 0.021 
Callibaetis 1 37.8 18.1 8.26 0.029 
Macrothemis 1 34.8 17.8 7.52 0.029 
Paracloeodes 1 33 16.1 8.11 0.047 
Meridiallaris/Massartellopsis 3 71.4 18.8 8.77 0.001 
Trepobates 3 44.5 27.4 9.05 0.05 
Cylloepus 3 40.9 19.8 8.22 0.03 
Athericidae 3 37.9 14.6 7.16 0.015 























Table IV. Spearman rank correlations of biotic metrics with classes of disturbance level. Only metrics 
with a statistically significant value (p < 0.05) are shown, and the highest Spearman coefficients are 
highlighted in bold. 


Disturbance level p 
Shannon Diversity -0.458 0.010 
Simpson Diversity -0.37 0.043 
% EPT Individuals 0.394 0.031 
Ephemeroptera Ind 0.454 0.011 
% Ephemeroptera Ind 0.488 0.006 
% Ephemeroptera Таха 0.403 0.027 
Elmidae Ind -0.532 0.002 
Elmidae Taxa -0.487 0.006 
% Elmidae Ind -0.552 0.001 
% Elmidae Taxa -0.433 0.016 
Trichoptera Taxa -0.532 0.002 
% Trichoptera Таха -0.372 0.043 
N Trichoptera -0.519 0.003 
Plecoptera Ind -0.629 0.000 
% Plecoptera Ind -0.683 0.000 
Plecoptera Taxa -0.574 0.000 
% Plecoptera Taxa -0.448 0.013 
EIPT Ind -0.42 0.020 
EIPT Taxa -0.635 0.000 
% EIPT Ind -0.394 0.031 
% EIPT Taxa -0.561 0.001 
Coleoptera Ind -0.474 0.008 
% Coleoptera Ind -0.474 0.001 


Table V. Indicator taxa for disturbance level (WP: well preserved, М=11; Di: disturbed, N=22) obtained 
by combining abundance and frequency data according to Dufréne & Legendre's (1997) method. Only 
taxa with a statistically significant indicator value (larger than expected by chance, p « 0.05) are shown. 
Indicator values range from 0 to 100 (perfect indicator). 









































Taxa Disturbance level Indicator value Mean 50 р 
Anacroneuria МУР 74.8 38.6 7.9 0.001 
Ambrysus WP 64.8 38.1 8.81 0.009 
Calamoceratidae WP 55.5 37.1 8.3 0.036 
Baetodes Di 51 29.1 7:9 0.021 
Corydalus МУР 49.5 327 7.08 0.038 
Teinopodagrion meridionale WP 42.1 9.1 7.11 0.01 
Cylloepus WP 41.3 8.6 6.82 0.01 
Brechmorhoga nubecula МУР 40.9 8.8 6.99 0.014 
Cernotina WP 39.8 8.6 6.8 0.02 
Marilia WP 34.8 8.7 6.95 0.023 
Helichus WP 31.5 7.1 6.67 0.042 
Gymnochtebius МУР 27.3 0.9 4.9 0.035 
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Table VI. Indicator taxa for factors of habitat disturbance (A: absence of, P: presence of) obtained by 
combining abundance and frequency data according to Dufréne & Legendre's (1997) method. Only 
taxa with a statistically significant value (larger than expected by chance, p « 0.05) are shown. Indicator 


values range from 0 to 100 (perfect indicator). 
















































































Taxa Variable Indicator value Mean SD p 

Progomphus P. channelization (N22) 89.7 29.1 12.47 0.011 
Thraulodes P. channelization (N22) 82.2 52.1 14.82 0.028 
Baetodes P. crops (N=5) 60 32.4 10.46 0.025 
Berosus 9 P. crops (N=5) 40 9.9 5.45 0.015 
Tropisternus P. crops (М=5) 38.6 13.1 6.7 0.046 
Baetodes P. roads (М=12) 56.7 28.7 7.58 0.003 
Baetodes P. cattle (N=13) 49.6 28.6 737 0.015 
Thraulodes P. logging (М=15) 60 39.4 7.56 0.017 
Baetodes P. logging (N-15) 47.1 28.4 7.23 0.023 
Phyllocycla P. logging (М=15) 26.7 3:7 5.6 0.046 
Anacroneuria A. roads (N=21) 60.2 38 7.69 0.008 
Teinopodagrion meridionale A. roads (N=21) 33.3 8.4 6.5 0.045 
Nectopsyche A. roads (N=21) 33.3 8.8 6.76 0.05 
Anacroneuria A. cattle (N=20) 59 38.1 7.48 0.016 
Nectopsyche A. cattle (N=20) 35 8.1 6.76 0.035 
Anacroneuria A. logging (N=18) 57.1 37.9 7.05 0.015 
Laccophilus A. logging (N=18) 41.4 21.7 6.74 0.015 
Teinopodagrion A. logging (N=18) 31.5 7.9 6.32 0.046 
Marilia A. logging (N=18) 30.6 8.3 6.36 0.038 
Americabaetis A. garbage (N=27) 58.8 39.7 9.9 0.024 
Callibaetis A. garbage (N=27) 41.8 20.7 8.34 0.027 
Grumichella A. garbage (N=27) 29.5 67 7.38 0.043 




















analysis overlaid onto ordination diagrams 
can be characterized according to the same 
variables: groups 4 and 5 including the 
totality of foothill rainforest localities, of lower 
elevation, higher water temperature and 
wider channel (Figs. 3.1-3.2), are separated 
by percentages of undercut banks, large 
woody debris in the channel and coarse 
gravel in the substrate (higher for group 4; 
Fig. 3.2). Groups 1 to 3, including all 
mountain forest localities of higher elevation 
with lower water temperature and narrower 
channels (Figs. 3.1-3.2), are ordered on the 
basis of latitude and percentage of cobble and 
large woody debris (Figs. 3.1-3.3). 
Disturbance level was not significantly 
correlated with ordination axes, and neither 
was pH in the second NMS analysis (run 
excluding 8 localities lacking pH data and 
taxa present only in them). 

Indicator analysis for altitudinal floors 


yielded ten taxa with a significant indicator 
value (Table 11), of which the best are 
Americabaetis (Ephemeroptera, Baetidae) for 
foothill rainforest streams, and Meridiallaris/ 
Massartellopsis(Ephemeroptera, Leptophlebiidae) 
for mountain forest streams. In accordance 
with the results shown by the ordination 
analyses, no indicator taxa were found for 
the intermediate altitudinal level (mountain 
rainforest). 


Analysis of disturbance 


Disturbance gradient was not significantly 
correlated with any of the natural 
environmental variables measured. 

Comparison of sites in well preserved 
areas (11 localities; mean inner distance of 
0.799) with sites in disturbed areas (22 
localities; mean inner distance of 0.852) by 
means ofthe MRPP test showed the difference 
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in their aquatic insect assemblage 
composition to be significant (observed delta: 
0.835; expected delta: 0.843; variance of 
delta: 0.17%; T: -2.088; p: 0.035; A: 0.011). 

Richness, diversity and several biological 
metrics (ѓе. Plecoptera individuals, 
proportion of Elmidae individuals, 
Trichoptera taxa, proportion of 
Ephemeroptera individuals and EIPT taxa, Fig. 
4) were found to be significantly different 
among sites in well preserved and in 
disturbed areas, not so EPT taxa (Fig. 4). 

Both Shannon and Simpson diversity 
indexes were significantly correlated with the 
disturbance gradient. Ephemeroptera metrics 
were either directly correlated with it 
(number of individuals and proportions of 
individuals and of taxa), increasing with 
higher disturbance levels, or not correlated 
(number of taxa), and the traditional EPT taxa 
index was accordingly not sensitive to this 
gradient (Table IV). Correlations of metrics 
calculated for Plecoptera, Elmidae, 
Trichoptera, and combination of these three 
taxa (EIPT) were highly significant. Among 
these, the best correlated with the disturbance 
gradient were proportion and number of 
Plecoptera individuals (r -0.68 and -0.63), 
and EIPT taxa (r -0.635), all of them beyond 
the 0.001 significance level (Table IV). 

Twelve taxa with a significant indicator 
value for disturbance level were identified 
(Table V), with Anacroneuria (Plecoptera, 
Perlidae) as the best indicator for sites in well 
preserved areas. Both positive and negative 
indicators significantly correlated with 
individual disturbance factors were found 
(Table VI). 


DISCUSSION 


Of the 55 families of aquatic insects 
found, 51 correspond to benthic insects, and 
the number of families per site is similar to 
that found in other studies of benthic 
macroinvertebrates in Andean mountain 
streams (Salinas et a/., 1999 in Bolivia, 
Fernández, et al., 2002 in Tucumán, and 
Miserendino, 1995 in Patagonia). Diversity 
values are also comparable to those obtained 
by Fernández et al., 2001 in a river in the 
Yungas of NW Argentina. The highest richness 


and diversity values observed in protected 
areas emphasize the importance of habitat 
conservation, and show that these areas are 
fulfilling their role in preserving the aquatic 
insect biodiversity of the Yungas. 

Altitude was the strongest variable shaping 
studied communities of aquatic insects, but 
in accordance to the results of Maldonado & 
Goitía (2003) for benthic macroinvertebrates 
in Yungas streams of Bolivia, communities 
were not perfectly correlated with the three 
altitudinal floors; those from the intermediate 
altitude level (mountain rainforest) were not 
separated as a distinct group but were 
interspersed between low and high altitude 
level assemblages. The significant correlation 
of channel width and substrate size with 
community composition observed here was 
found also in studies of benthic 
macroinvertebrates of lotic Yungas 
environments (Fernández et a/., 2001; Moya 
et al., 2003; Oller & Goitía, 2005), and of 
Patagonian rivers (Miserendino & Pizzolón, 
2000; 2003). 

An understanding of relationships 
between anthropogenic land use and the 
ecological status of streams can be 
complicated by covariation between 
anthropogenic and natural gradients (Allan, 
2004). In this study no covariation between 
both gradients was found, and disturbance 
level was not significantly correlated with any 
ordination axis, thus indicating that natural 
variables are the primary factors influencing 
composition of aquatic insect communities 
of Yungas streams. 

In a study of four sites at the Lules River 
in the Yungas of Tucumán province with 
values of physical and chemical parameters 
and of biotic indexes (BMWP, ASPT, EPT) 
within the ranges of non-polluted water, 
Fernández et al. (2002) found that neither 
BMWP nor EPT indexes reflected the riparian 
disturbance observed in two of these sites. 
Approximations to BMWP and EPT in the 
present study did not correlate with riparian 
disturbance either, but metrics comprising 
Elmidae, Trichoptera and Plecoptera 
presented highly significant correlations with 
it (Table IV). Although Coleoptera and 
Elmidae are not always included in biological 
indexes (Lenat, 1993; Prat et al., 1999; Klemm 
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et al., 2003), their potential as bioindicators 
and as part of biological metrics has been 
studied in Europe (Richoux & Forestier, 1989; 
García-Criado & Fernández-Alaez, 1995; 
2001) and the USA (Ode et a/., 2005). Fossati 
et al. (2001) found that Elmidae were among 
the most affected insects by sediment releases 
due to road construction in eight sites along 
an Andean river in Bolivia, and this study 
suggests that they would be a useful 
component in an index sensitive to the 
disturbance level at the local scale for Yungas 
streams. 

EPT metrics have been found to be 
negatively correlated with disturbance in 
several studies (/.e. Sponseller et a/., 2001; 
Klemm et a/., 1990; 2003; Miserendino, 
2004; Morais et al., 2004). The positive 
correlation of Ephemeroptera metrics with the 
disturbance gradient found here accounts for 
the lack of correlation of EPT with that 
gradient. A number of environmental stressors 
can have positive influence at low to 
moderate concentrations; riparian thinning 
and low intensity agriculture increase light, 
nutrients, and water temperatures, leading to 
an increase of periphyton biomass (Allan, 
2004), which can benefit some species in 
spite of their negative effects (/.е. decrease of 
bank stability, inputs of litter and wood, and 
retention of nutrients, and increase of bank 
and channel erosion and input of pesticides; 
Gregory et al., 1991). Some Ephemeroptera 
taxa benefited from disturbance, and this 
explains the positive correlation of 
Ephemeroptera metrics with the disturbance 
gradient; Baetodes (Ephemeroptera, Baetidae) 
was found with significantly higher 
abundance and frequency in disturbed sites 
(Table V), and was significantly associated to 
the presence of particular disturbance factors, 
as was Thraulodes (Ephemeroptera, 
Leptophlebiidae) (Table VI). 

The high values found for Progomphus 
(Odonata, Gomphidae) and Thraulodes 
(Ephemeroptera, Leptophlebiidae) as 
indicators of channelization could be 
accounted for by changes in the 
representation of their preferred habitats 
resulting from both direct and indirect effects 
of channelization. Many channel or stream 
bank stabilization structures provide 


increased instream habitat for certain aquatic 
species; Sandheinrich & Atchison (1986) 
reported increases in densities of epibenthic 
insects within revetments and stone dike areas 
and more suitable substrate for bottom- 
dwelling insects in revetment areas, which 
could explain the high value obtained for 
Thraulodes. One of the secondary instream 
habitat alteration effects of channelization is 
the trapping of large quantities of sediment 
(EPA, 2006), which could favor an increase 
in populations of Progomphus larvae 
(burrowers in sediments). In spite of the high 
indicator and significance values that the 
analysis shielded for these two taxa, further 
study is required to verify their validity 
because they were based on only two sites 
showing evidence of channelization. 

Results suggest that composite measures, 
such as a multimetric index including 
richness, diversity and biological metrics such 
as an EIPT (Elmidae, Plecoptera, Trichoptera) 
index, could be useful in detecting overall 
disturbance at a local scale in Yungas streams. 
Further studies are needed to develop and 
test such an index. In order to carry out 
sound management and restoration actions 
effectively, it is necessary to diagnose causes 
as well as to assess disturbance; because of 
the aggregated nature of multimetric indexes, 
they might be less easily interpreted than the 
behavior of individual variables (Watzin & 
McIntosh, 1999). Although indicator values 
for individual disturbance factors were in 
general low in this study, these taxa are worth 
investigating in future studies analyzing their 
tolerance values to specific stressors resulting 
from disturbance factors (Black et a/., 2004; 
Yuan & Norton, 2004), thus providing the 
diagnostic tools necessary to identify and 
monitor disturbance factors acting upon 
Yungas streams. 
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